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Abstract 

' We explore the possibility of extrapolating state of the art lattice QCD calculations of nucleon 

form factors to the physical regime. We find that the lattice results can be reproduced using the 
— . ■ Light Front Cloudy Bag Model by letting its parameters be analytic functions of the quark mass. 

We then use the model to extend the lattice calculations to large values of Q 2 of interest to current 
[ and planned experiments. These functions are also used to define extrapolations to the physical 

value of the pion mass, thereby allowing us to study how the predicted zero in Ge(Q 2 )/Gm(Q 2 ) 

varies as a function of quark mass. 
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I. INTRODUCTION 



The electromagnetic form factors of the nucleon are an invaluable source of information 
on its structure For example, observing their fall as Q 2 increases from zero revealed the 
finite extent of the nucleon, and measuring the Sachs electric form factor of the neutron, 
G% that it has a positive core surrounded by a long-range, negative tail IH IH E[ . In the 

last few years particular interest has focused on the ratio of the electric and magnetic form 
factors of the proton, Ge/Gm, where recoil polarization data 0, S| have revealed a dramatic 
decrease with Q 2 - in contrast with earlier work based on the Rosenbluth separation. These 
data have allowed one of us to deduce a fascinating spin dependence of the shape of the 
nucleon 0. 

While the behavior of Ge/Gm with Q 2 was anticipated in some models (e.g. see Ref. (Io| . 
[ill]), there is no consensus as to which explanation best represents how QCD works. Direct 
guidance from QCD itself would be most valuable and for that purpose lattice QCD rep- 
resents the one and only technique by which one can obtain non-perturbative solutions to 
QCD. 

The QCDSF Collaboration recently presented lattice QCD simulations for the form fac- 
tors of the nucleon over a wide range of values of momentum transfer l3| ■ While these were 



based on the quenched approximation, with an unsophisticated action, several lattice spac- 
ings were chosen with the smallest being around 0.05 fm (/3 = 6.4) and at present these are 
the state of the art. The quark masses used in the simulations correspond to pion masses in 
the range (0.6, 1.2) GeV. Therefore one needs to parametrize the form factors as a function 
of pion mass and extrapolate to the physical value before comparing these lattice results 
with the experimental data. 

At Q 2 = there have been a number of studies of the chiral extrapolation of baryon 
magnetic moments 0, 3 2 3- However, there is no model independent way to respect 



the constraints of chiral symmetry over the range of Q 2 and m n required by the QCDSF 
data. Instead, at finite Q 2 , one has been led to study various phenomeno logical parameteri- 
zations |rH |. which have at least ensured the correct leading order non-analytic structure as 
m n — > 0. Our purpose here is three- fold. First, we wish to use the lattice data to investigate 
whether a particular quark model is capable of describing the properties of the nucleon in 
this additional dimension of varying - an important test which any respectable quark 
model should satisfy 1 . Second, having confirmed that the model is consistent with the lattice 
data over the range of rn n noted earlier, we use the model to extrapolate to large values of 
Q 2 (for lattice values of m^). Third, we also use the model to extrapolate to the physical 
pion mass. 

The model which we consider here is the light front cloudy bag model (LFCBM) [l9| . 
which was developed as a means of preserving the successes of the original cloudy bag model 
0], while ensuring covariance in order to deal unambiguously with modern high energy 
experiments. The light front constituent quark model, upon which it is built ^(j , predicted 
the rapid decrease of Ge/Gm with Q 2 and, as the pion cloud is expected to be relatively 
unimportant at large Q 2 , this success carries over to the LFCBM Furthermore the 

LFCBM corresponds to a Lagrangian built upon chiral symmetry, so it can be extended to 



1 Just as the study of QCD as a function of N c has proven extremely valuable, so the study of hadron 
properties as a function of quark mass, using the results of lattice QCD calculations, undoubtedly offers 
significant insight into QCD, as well as new ways to model it 18]. 
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the limit of low quark mass as well as low and high Q 2 . 

The outline of the paper follows. In Sect. 2 we briefly review the LFCBM. In Sect. 3 we 
present the lattice QCD data, explain the fitting procedure and present the results. Sect. 4 
contains some concluding remarks. 

II. REVIEW OF THE LFCBM 

The light front cloudy bag model (LFCBM) respects chiral symmetry and Lorentz invari- 
ance and reproduces the four nucleon electromagnetic form factors. Therefore it is reasonable 
to try to use it to extrapolate the form factors computed using lattice QCD to the physical 
pion mass. We begin by briefly introducing the key features of the LFCBM. 

The LFCBM is a relativistic constituent quark model incorporating the effect of pion- 
loops, key features motivated by chiral symmetry. The light-front dynamics is employed to 
maintain the Poincare invariance, and one pion-loop corrections are added to incorporate 
significant pion cloud effects (particularly in the neutron electric form factor and magnetic 
moments) as well as the leading non-analytic behavior imposed by chiral symmetry. In 
light-front dynamics the fields are quantized at a fixed "time"= r = x° + x 3 = x + . The light 
front time or r-development operator is then P° — P 3 = P . The canonical spatial variable 
is x~ = x° — x 3 , with a canonical momentum P + = P° + P 3 . The other coordinates are 
xi and Pj_. The relation between the energy and momentum of a free particle is given by 
p~ = (p 2 ± + m 2 )/p + , with the quadratic form allowing the separation of center of mass and 
relative coordinates. The resulting wave functions are frame independent. The light front 
technique is particularly relevant for calculating form factors because one uses boosts that 
are independent of interactions. 

Our goal is to calculate the Dirac Pi and Pauli P2 form factors given by: 

(N, \'p'\J»\N,\p)=u x ,(p') 

The momentum transfer is g M = ( p' — p) M , Q 2 = —q 2 and P* is taken to be the electro- 
magnetic current operator for a free quark. For Q 2 = the form factors Pi and F 2 are, 
respectively, equal to the charge and the anomalous magnetic moment k in units of e and 
e/{2Mjsi) , and the magnetic moment is \i — Pi(0) + F 2 (0) — 1 + k. The evaluation of the 
form factors is simplified by using the so-called Drell-Yan reference frame in which q + = 0, 
so that Q 2 = q\ = q\. If light-front spinors for the nucleons are used, the form factors can 
be expressed in terms of matrix elements of the plus component of the current [2jJ: 

Pi(Q 2 ) = (JV,t \ J + \N^}, and QF 2 (Q 2 ) = (-2M N )(N,] |J + |iV,|). (2) 

The form factors are calculated using the "good" component of the current, J + , to suppress 
the effects of quark-pair terms. Finally, we note that in our fits we will use the Sachs form 
factors, which are defined as 

Q 2 

Ge = Pi — , , j9 F2, Gm = Pi + P2 (3) 

The next step is to construct the bare (pionless) nucleon wave function which is a 
symmetric function of the quark momenta, independent of reference frame, and an eigenstate 
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of the canonical spin operator. The commonly used ansatz is: 



*(pi) = $(M 2 )m( Pi )m(p 2 )«(p3)^Pi,P2,P3), (4) 

Pi Pi^i; Tj) 

where ^ is a spin-isospin color amplitude factor, the pi are expressed in terms of relative 
coordinates, the u(pi) are Dirac spinors and $ is a momentum distribution wave function. 
The specific form of ip is given in Eq. (12) of Ref. and earlier in Ref. [llj]. This 
is a relativistic version of the familiar SU(6) wave function, with no configuration mixing 
included. The notation is that p.; = (pf, Pij_). The total momentum is P = pi + p 2 + P3, 
the relative coordinates are £ = Px/(pt +pt), V = (pf + pt)/P + 1 an d kj_ = (1 — £)pi_i_ — 
£p 2 j_, Kj_ = (1 — 77)(pij_ + P2±) — ?7P3_l- In computing a form factor, we take quark 3 to be 
the one struck by the photon. The value of 1 —rj is not changed (q + = 0), so only one relative 
momentum, Ki is changed: K' ± = Kx — ?7q_i_. The form of the momentum distribution 
wave function is taken from Schlumpf [22]: 



$(M ) = iV(M 2 + /3 2 ) 7 , (5) 
with Mq the mass-squared operator for a non-interacting system: 

2 K\ k\ + M 2 M 2 

Schlumpf's parameters were (3 = 0.607 GeV, 7 = —3.5, M = 0.267 GeV, where the 
value of 7 was chosen so that Q 4 Gm(Q 2 ) is approximately constant for Q 2 > 4 GeV 2 , in 
accord with experimental data. The parameter (3 helps govern the values of the transverse 
momenta allowed by the wave function $ and is closely related to the rms charge radius. 
The constituent quark mass, M, was primarily determined by the magnetic moment of the 
proton. We shall use different values when including the pion cloud and fitting lattice data. 

A physical nucleon can sometimes undergo a quantum fluctuation so that it consists of 
a bare nucleon and a virtual pion. In this case, an incident photon can interact electro- 
magnetically with a bare nucleon, Fig. [TJi, with a nucleon while a pion is present, Fig. ^d, 
or with a charged pion in flight, Fig. ^p. These effects are especially pronounced for the 
neutron Gb[6], at small values of Q 2 . The tail of the negatively charged pion distribution 
extends far out into space, causing the mean square charge radius, i? 2 , to be negative. The 
effects of the pion cloud need to be computed relativistically if one is to confront data taken 
at large Q 2 . This involves evaluating the Feynman diagrams of Fig. fusing photon-bare- 
nucleon form factors from the relativistic model, and using a relativistic 7r-nucleon form 
factor. The resulting model is defined as the light-front cloudy bag model LFCBM [19]. 
The light-front treatment is implemented by evaluating the integral over the virtual pion 
four- momentum , kj_, by first performing the integral over k~ analytically, re-expressing 
the remaining integrals in terms of relative variables (a = k + /p + ), and shifting the rel- 
ative x variable to to simplify the numerators. Thus the Feynman graphs, Fig. 
are represented by a single r-ordered diagram. The use of J + and the Yan identity [2c 
Sf(p) — u {Pi s )u(p, s)/{p 2 — m 2 -Me) +7 + /2p + allows one to see that the nucleon current 
operators appearing in Fig^ act between on-mass-shell spinors. 

The results can be stated as 

Fia(Q 2 ) = Z \f£\q 2 ) + F tbQ (Q 2 ) + F lCQ (Q 2 )] , (7) 



4 



N I N N , N^N N N f N N 

N 

k 

(a) (b) (c) 

FIG. 1: Diagrams 



where i — 1,2 denotes the Dirac and Pauli form factors, a = n,p determines the identity 
of the nucleon, and F^(Q 2 ) are the form factors computed in the absence of pionic effects. 
The wave function renormalization constant, Z, is determined from the condition that the 
charge of the proton be unity: Fi p (Q 2 = 0) = 1. For illustration we start with the calculation 
of the neutron form factors. Then, evaluating the graph in Fig. gives 

F lbn (Q 2 ) = g 2 j\aa j ^fl^W \ a)R N (L^ \ a) 

[(fS\Q 2 ) + F^(Q 2 )/2)(a 2 (M 2 - Q 2 /4) + L 2 ) - (f£\Q 2 ) + F^(Q 2 ) /2))(a 2 Q 2 /2)], (8) 

pX p i2 j 

F 2bn {Q 2 ) = -9o J daa J —-R N (L^ \ a)R N (L^ \ a) 

K ) (Q 2 ) + ^ ) (Q 2 ))(2« 2 M^ 



2 

(F£\Q 2 ) + ^\Q 2 )))(a 2 M 2 (l - Q 2 /AM 2 ) + (L 2 X - L 2 y ))] (9) 



where go is the bare 7rN coupling constant, and the renormalized coupling constant Zg^ = 
g 2 /47T = 13.5 ,Lf } = L± ± aq±/2, a = ^, D N (k 2 ± ,a) = M 2 a 2 + k\ + /i 2 (l - a), and 

R N (k 2 ± , a) = F J^i'ay Th e vrN form factor is taken as @ 

F nN (ki, a) = e -(°N(kl,a)/2a(l-a)Al)^ (1Q) 



and maintains charge conservation [271] . The constant A is a free parameter, but very large 
values are excluded by the small flavor asymmetry of the nucleon sea. 

From Eqns. (JHJ) and © we see that each term in the nucleon current operator contributes 
to both Fi and F 2 . The evaluation of graph ^ yields 



F lcn (Q 2 ) = -9 2 F n (Q 2 ) J 1 daa j 0-R(K™ \ a)R(K^ 2 , 



K i + M 2 a 2 -(l-a) 2 ^- 

4 

fill' 



F2UQ 2 ) = -g 2 (2M 2 )F w (Q 2 ) jf * daa 2 {\ - a) j |^(K (+) a)R(K^ 2 , a) (12) 
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where Kf } = Kj_ ± (1 - a)q ± /2 2 

The proton form factors can be obtained by simply making the replacements n — > p in 
Eqs. (I8I9|) and — g 2 — > g 2 . in Eqs. (jllll2J) . The change in sign accounts for the feature that 
the tt~ cloud of the neutron becomes a tt + cloud for the proton. The mean-square isovector 
radii (r 2 )J , computed using Eqs. Q, and then taken to the chiral limit at low-CP, have the 
same singular log terms as those of the relativistic results of Beg and Zepeda |26j ]. 

The LFCBM was defined by choosing four free parameters: m, (3, 7, A so as to best 
reproduce the four experimentally measured electromagnetic form factors of the nucleon 
Til l . In the present work, the most relevant of these parameters will be varied to reproduce 



lattice data, and the resulting dependence on the quark mass and lattice spacing used to 
extrapolate to the physical region. 



III. FITTING THE QCDSF FORM FACTORS AND EXTRAPOLATING TO THE 
PHYSICAL PION MASS. 

In this section we discuss the fitting procedure used to parametrize the nucleon form 
factors calculated in lattice QCD. We use data produced by the QCDSF Collaboration 
T2l | and employ the LFCBM to calculate the corresponding form factors, varying the model 



parameters to find the best-fit to the different sets of lattice data obtained for different values 
of the current quark mass, m q . The behavior of the fitting parameters is then represented by 
a polynomial function of the quark mass m q . This polynomial fit in m g , or equivalently in 
pion mass squared, m 2 , can then be used to extrapolate the values of the fitting parameters to 
the physical pion mass. Nucleon form factors for the physical pion mass are then calculated 
using the extrapolated values for the model parameters. In the following few subsections a 
more elaborate explanation is given and the results are presented. In section lHI Al we describe 
the available data and the analysis procedure used to extract the quantities necessary for 
further fits. In sec. IIIIBI we describe the details of the fitting and extrapolation process and 
in the sec. IHI CI we present the nucleon form factors resulting from the extrapolation to the 
physical pion mass and make comparisons with experiment. 



A. QCDSF Data and Its Analysis 

The form factor calculations in Ref.[l2^ were carried out for three different values of 
the lattice spacing, a = {0.47, 0.34, 0.26} GeV -1 . For each value of a several sets of pion 
(or equivalently nucleon) masses were considered. For each mass set Dirac and Pauli form 
factors for both the proton and neutron were calculated at several values of Q 2 . The typical 
range for the pion mass used varied from 1.2 GeV to 0.6 GeV, with the corresponding nucleon 
mass ranging from approximately 2 GeV to 1.5 GeV. The typical range for Q 2 was 0.6 GeV 2 
to 2.3 GeV 2 . 

The LFCBM is basically a relativistic constituent quark model, so we need to relate the 
model constituent mass of Eq. © to the masses of the nucleon and pion. To do so we use 



2 These formulae are slightly different from those of Ref. [19j . This leads to slight changes in the parameters 
that will be discussed elsewhere. 
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the approach of Ref. 0, ( Eq. (8)) 



cm 



phys 

1 ™2 



M = M X + ^f^X, (13) 



[m 



where M x is the constituent quark mass in the chiral limit, m v ^ ys is the current quark mass 
and c is of order 1. In the study of octet magnetic moments in the AccessQM model of 
Ref. 0, the best fit value for M x was 0.42 GeV, while for cm p q hys it was 0.0059 GeV. 



B. Lattice Data Fit and Extrapolation 

The first step in our extrapolation of the lattice results to the physical quark mass is to 
fit the lattice results for each quark mass m q by adjusting the parameters of the LFCBM 
calculation. For that purpose two fitting parameters were chosen. The first parameter is M x 
in Eq. (j!3|) . which determines the constituent quark mass. This parameter was varied for 
each lattice spacing separately, since some dependence upon lattice spacing was anticipated. 
The second parameter is the internal parameter, 7, in the nucleon wave function Eq. (j3J), 
which is varied separately for each pion (or equivalently nucleon) mass. For convenience, we 
express all magnetic form factors Gm in "physical" units of e/2M^ hysical . Since the LFCBM 
uses the mass of the p-meson included in the pion electromagnetic form factor, we need the 
extrapolated value for its mass. We use the simple fitting function from Ref. |28| : 

m p = c + cxml, (14) 

with c = 0.776 GeV and c x = 0.427 GeV" 1 . 

A function representing the x 2 f° r the deviation between the lattice data and the values 
calculated using the LFCBM was constructed and minimized by varying the fitting param- 
eters. Changing the value of M x causes the the calculated form factors to move up or down 
by an amount approximately independent of Q 2 , thereby causing a relatively small change 
in x 2 ■ Therefore a simple grid variation for that parameter was employed, with grid bound- 
aries M x G [0.15, 0.45] GeV, and step size of 5m x = 0.01 GeV. As for the parameter, 7, 
the variation of x 2 was much stronger and the Minuit package of CERN's Root framework 
[29l | was used for the minimization. At first the boundaries for 7 were set to keep it in 
the physical region, but successful boundless runs were also performed in order to confirm 
the true minimum and error sizes. The pion masses used in the lattice calculation are very 
large, and the resulting pionic effects are very small. Therefore the value of A could not be 
determined from lattice data and its value was held fixed at A = 0.58 GeV Similarly, varying 
(3 did not change the description of the lattice data, so it was held fixed at (3 = 0.607 GeV /c. 
The resulting fits are in good agreement with data, as one can see in Figs. El HI and [SJ 
The best-fit values of the parameters are shown in Table I. The figures show results for the 
smallest lattice spacing, a = 0.26 GeV -1 , but the reproduction of lattice data is equally 
successful for larger values of a. 
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. 2: (Color online) LFCBM fit to QCDSF data for G^(in units of e) for a lattice spacing 
0.26 GeV" 1 , M P = 1.80 GeV and = 0.93 GeV 
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FIG. 3: (Color online) LFCBM fit to QCDSF data for G p M (in units of e/(2M^ hysical )) for a lattice 
spacing a = 0.26 GeV" 1 , M P = 1.80 GeV and m w = 0.93 GeV 
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FIG. 4: (Color online) LFCBM fit to QCDSF data for G^(in units of e) for a lattice spacing 
a = 0.26 GeV" 1 , M P = 1.80 GeV and = 0.93 GeV 
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FIG. 5: (Color online) LFCBM fit to QCDSF data for G%(in units of e/(2M^ hysical )) for a lattice 
spacing a = 0.26 GeV" 1 , M P = 1.80 GeV and m n = 0.93 GeV 
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TABLE I. Lattice data and LFCBM fitting parameters. (All expressed in powers of GeV.) 



a 


m n 


M N 


M x 


7 


0.47 


1.146 


2.062 


0.390(5) 


-6.12(7) 


0.47 


1.068 


1.981 


0.390(5) 


-5.67(6) 


0.47 


0.873 


1.746 


0.390(5) 


-4.95(9) 


0.47 


0.752 


1.567 


0.390(5) 


-4.78(12) 


0.47 


0.638 


1.503 


0.390(5) 


-4.67(15) 


0.47 


0.135 


0.938 


0.390(5) 


-4-79(46) 


0.34 


1.201 


2.141 


0.280(5) 


-5.03(6) 


0.34 


1.035 


1.933 


0.280(5) 


-4.37(5) 


0.34 


0.881 


1.732 


0.280(5) 


-4.99(5) 


0.34 


0.706 


1.522 


0.280(5) 


-3.51(6) 


0.34 


0.135 


0.938 


0.280(5) 


-2.91(29) 


0.26 


1.237 


2.202 


0.210(5) 


-4.78(5) 


0.26 


1.092 


2.028 


0.210(5) 


-4.14(7) 


0.26 


0.925 


1.802 


0.210(5) 


-3.69(5) 


0.26 


0.744 


1.600 


0.210(5) 


-3.09(6) 


0.26 


0.580 


1.379 


0.210(5) 


-3.01(13) 


0.26 


0.135 


0.938 


0.210(5) 


-2.41(22) 



The next step is to extrapolate the fitting parameters to the physical quark mass. This is 
done using the assumption that the parameters vary smoothly as functions of the quark mass, 
and the fact that m q ~ m 2 over the mass range investigated. We limited the extrapolation 
function to a low order polynomial in m 2 . The resulting fits for two lattice spacings are 
presented in Figs. |H1 and from which we see that the fitting function provides a very 
accurate representation of the values obtained from lattice data. The fitted values of 7 and 
the extrapolation to the physical value of m n , with their corresponding errors, are shown in 
Figs. Eland[7| 

In our calculations, M x has a very weak dependence on the pion mass, but it has a rather 
strong dependence upon the lattice spacing. As we see in Table I and Figs. I2E very good 
fits to the lattice data are obtained even without varying M x for each quark mass. By 
contrast, Fig. |S] and Table I show rather dramatic variation of M x for different values of 
the lattice spacing, a. This suggests that the larger values of the lattice spacing are rather 
far from the continuum limit and (at best) only the results for the smallest lattice spacing 
should be compared with experimental data. It would clearly be desirable to have new data 
at even smaller a, or using an improved action, known to provide a good approximation to 
the continuum limit. 

Use of the values of 7, M determined by the lattice data in the LFCBM defines a lattice 
version of the LFCBM. We may use this new model to compute the form factors at arbitrarily 
large values of Q 2 , thereby extending the kinematic range of the lattice calculations. The 
results are shown in Figs. El El ini an d El in Figs. El and El we show the corresponding 
plots of HqGe/Gm- 
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FIG. 6: (Color online) Polynomial extrapolation of 7 vs. m n for lattice spacing a = 0.47 GeV 
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FIG. 7: (Color online) Polynomial extrapolation of 7 vs. for lattice spacing a = 0.26 GeV 
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FIG. 8: (Color online) Variation of M x with lattice spacing a. The best-fit values of M x for physical 
m n using LFCBM and AcessQM models are presented as well. 



C. Results at the Physical Pion Mass and Comparison With Experiment 

We use the extrapolated values of 7 and M (Figs. EHHJ) to calculate the nucleon electric 
and magnetic form factors using the physical pion and nucleon masses. The resulting plots 
for Ge-i Gm and their ratios vs. Q 2 for both proton and neutron are shown in Figs. IToTOIH 
Figure El shows that our results are in more or less good agreement with the experimental 
data in the low-Q 2 region, but yield a slightly lower value of Q 2 for the zero cross-over point 
than that extrapolated from experiment j^Tj. A new analysis that includes an estimate of 
all of the effects of two photon exchange yields a zero-crossing value that is somewhat closer 
to ours [31] but future data will resolve this unambiguously. 

An alternative method of determining the value of the Q 2 for which Ge/Gm passes 
through zero at the physical pion mass is to fit the crossover values as a linear function of 
m 2 and extrapolate again to the physical pion mass. The resulting plot is shown in Fig. |^] 
This procedure yields approximately the same cross-over point as found in Fig. El 
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FIG. 9: (Color online) LFCBM calculations using parameters (Figs. H3[7|and|Hl) obtained by fitting 
the lattice results for the proton electric form factor, Ge, at lattice spacing a = 0.26 GeV^ 1 . 
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FIG. 10: (Color online) LFCBM calculations using parameters (Figs. El[I|and|HJ) obtained by fitting 
the lattice results for the proton magnetic form factor, Gm, at lattice spacing a = 0.26 GeV^ 1 . 
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FIG. 11: (Color online) LFCBM calculations using parameters (FigsElE|and|HJ) obtained by fitting 
the lattice results for the neutron electric form factor, Ge, at lattice spacing a = 0.26 GeV" 1 . 
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FIG. 12: (Color online) LFCBM calculations using parameters (FigsO30and|HJ) obtained by fitting 
the lattice results for the neutron magnetic form factor, Gm, at lattice spacing a = 0.26 GeV^ 1 . 
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FIG. 13: (Color online) LFCBM calculations using parameters (FigsJ5J [7| and |HJ) obtained by 
reproducing lattice results for the ratio of proton form factors, ^lqGe/Gm-, at lattice spacing 
a = 0.26 GeV- 1 . 
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FIG. 14: (Color online) LFCBM calculations using parameters (Figs|13 and IB} obtained by 
reproducing lattice results for the ratio of neutron form factors, ^oGe/Gm, at lattice spacing 
a = 0.26 GeV- 1 . 
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FIG. 15: (Color online) Extrapolated calculations for the proton electric form factor, Ge, for 
lattice spacing a = 0.26 GeV- 1 . The dashed and dotted curves show the upper and lower limits 
of variation of the calculated values due to the uncertainties of the parameters 7 and M x . 
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FIG. 16: (Color online) Extrapolated calculations for the proton magnetic form factor, Gm, for 
lattice spacing a = 0.26 GeV- 1 . The dashed and dotted curves show the upper and lower limits 
of variation of the calculated values due to the uncertainties of the parameters 7 and M x . 
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FIG. 17: (Color online) Extrapolated calculations for the ratio of proton form factors, [IqGe/Gmi 
for lattice spacing a = 0.26 GeV- 1 . The dashed and dotted curves show the upper and lower limits 
of variation of the calculated values due to the uncertainties of the parameters 7 and M x . 
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FIG. 18: (Color online) Extrapolated calculations for the neutron electric form factor, Ge, for 
lattice spacing a = 0.26 GeV- 1 . The dashed and dotted curves show the upper and lower limits 
of variation of the calculated values due to the uncertainties of the parameters 7 and M x . 
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FIG. 19: (Color online) Extrapolated calculations for the neutron magnetic form factor, Gm, for 
lattice spacing a = 0.26 GeV- 1 . The dashed and dotted curves show the upper and lower limits 
of variation of the calculated values due to the uncertainties of the parameters 7 and M x . 
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FIG. 20: (Color online) Extrapolated calculations for the ratio of neutron form factors, hqGe/Gm-, 
for lattice spacing a = 0.26 GeV' 1 . The dashed and dotted curves show the upper and lower limits 
of variation of the calculated values due to the uncertainties of the parameters 7 and M x . 
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FIG. 21: (Color online) Extrapolation of Qc ross where Ge/Gm passes through zero for lattice 
spacing a = 0.26 GeV^ 1 . 



IV. DISCUSSION 

Our study of the form factors calculated using the LFCBM with parameters determined 
by lattice data and by extrapolation to the physical pion masses yields very interesting 
results. The ratio G P E /G P M passes through zero for all of the calculations. The main variation 
of the position of the crossover between the fitting curves shown in Figs. EH and O comes 
from the variation of the nucleon mass, and not the variation of 7. Even though for the 
physical pion mass, the ratio varies rapidly as a function of 7 in the region 7 ~ —2, the 
function Ge/Gm for the neutron has a turning point at about 7 ~ —2.3. We shall explain 
these features using the LFCBM. 

Let us express the ratio Ge/Gm in terms of Pauli and Dirac form factors, F 1 and F 2 , 
respectively, using Eq. © 

G E _ Ft - Q 2 /{AMpF 2 1 + g 2 /(4M^) 

G M F, + F 2 l + F 1 /F 2 1 ' 

Consider first the values of Q 2 = {Qc ross ) where the ratio Ge/Gm for the proton passes 
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through zero for the set of calculations shown in Fig. ED Equation ()15|) tells us that 

Go™. = 4 <f (16) 

Now let us consider the formula for F ia (Q 2 ), Eq. (J7J). The second and third terms in 
Eq. (J7J) are only significant in the low — Q 2 region for physical pion masses. In the high — Q 2 
region, or for lattice calculations with high pion mass, these terms are vanishingly small. 
Indeed the numerical calculations support these statements, so we can neglect their contri- 
bution in the rest of the discussion. 

The corresponding formulas for and from Ref. are 

Fi°\Q 2 ) = j *Wo)*(M„) x fag" fa, p' 2 ) IxS^fri, P2)> (t Pal t Pa) (17) 



Q, ' rm 1 d2,hjk ([2Kj dV -&(M>MM ) x <xS ei (Pi,p' 2 )lxS ei (Pi,P2)) (T Psl I Pa) 



2M N j e(i-o^(i-^) 

(18) 

The $(M ) factors are wave functions of the form of Eq. (|5|). but using the lattice values 
of 7, M shown in FigsEHHI We stress that these two integrals differ only by the last factor, 
which gives the spin non-flip and spin-flip dependence of F^ and QF^ /2Mjy, respectively. 
At high Q 2 these matrix elements are each of order Q, causing the ratio QF 2 /Fi to be 
approximately constant. So we can express Q 2 Cross as 




Q 2 Cross = AM 2 N Uli (19) 



where denotes the integral for F^\ and L, denotes the integral for F^ . 

In the high-Q 2 region the ratio in Eq. (fT9*|l is approximately a constant, because the 
difference comes only from the overlap factors of the spin-dependent parts of the wave 
functions in the integrals (See fl9l]. j^). So the behavior of Qc ross is governed primarily by 
the factor M^. The linear variation of Q 2 Cross vs. presented in Fig. E21 shows the validity 
of this interpretation. 

We can also understand the behavior of Ge/Gm versus 7, by considering its role in the 
wave function. The factor 7 determines the size of the momenta appearing in the integrands 
of Eqs. ()17|) and ()18j) . The corresponding integrands differ by terms that are ratios of second 
order polynomials of the integration variables. For large absolute values of 7, the high 
momenta are cut off more strongly, so that the contribution of terms that cause differences 
between the integrals are not very significant. For small absolute values of 7 the integrals 
become more sensitive to those terms and we obtain a larger variation of the ratios of the 
integrals and hence the ratio Ge/Gm- 



V. CONCLUSION 

We have seen that the LFCBM can produce a very good description of the lattice QCD 
data for the nucleon form factors over a wide range of quark masses with a smooth, analytic 
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variation of the wave function parameter, 7, and the constituent quark mass, M. The pion 
cloud plays very little role in the mass range for which the lattice simulations have been 
made but it rapidly becomes more important as we approach the chiral limit. From the 
rather strong dependence of the form factors on the lattice spacing, a, it is not yet clear 
that we have obtained a good approximation to the continuum limit, but the form factors 
obtained at the smallest value of a are in reasonable agreement with experimental data in 
the low-Q 2 region for which the lattice simulations were made. 

At present the lattice simulations are limited to values of the momentum transfer at or 
below 2 GeV 2 and it is therefore a very big extrapolation to look at the behavior of the form 
factors in the region of greatest current interest. Nevertheless, the behavior of Ge/Gm which 
we find is particularly interesting. The ratio crosses zero for all values of the quark mass 
but the position where this happens varies over a very wide range of Q 2 . This variation can 
be understood almost entirely in terms of the variation of the corresponding nucleon mass, 
given that the ratio QFi/F 2 is approximately Q 2 independent in the model. We obtain the 
same value of Q 2 for the cross-over whether we extrapolate the position as a function of 
quark mass or simply evaluate the form factors at the physical pion mass, using the fitted 
dependence of the wave function parameters on pion mass. 

In the immediate future it is clearly very important to improve on the lattice data, both 
by ensuring that we really have a good approximation to the continuum limit (e.g., by using 



29 



a suitably improved action) and by extending the calculations to higher values of Q 2 . It 
would also be important to remove the need for quenching, even though that may not be 
such a limitation at large Q 2 . From the point of view of developing a deeper understanding 
of QCD itself it is important that the LFCBM is able to describe the present lattice data 
over such a wide range of masses. We would encourage a similar exercise for other models 
as a novel test of their validity. It remains to be seen whether the LFCBM has indeed been 
successful in predicting the behavior of the form factors at higher Q 2 and indeed whether it 
will match future experimental data. 
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